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Abstract 

The simplest extension of the MSSM that does not contradict LEP II experimental 
bound on the lightest Higgs boson mass at tan/3 ~ 1 is the modified Next-to- 
Minimal Supersymmetric Standard Model (MNSSM). We investigate the renor- 
malization of Yukawa couplings and soft SUSY breaking terms in this model. The 
possibility of 6-quark and r-lepton Yukawa coupling unification at the Grand Uni- 
fication scale Mx is studied. The particle spectrum is analysed in the vicinity of 
the quasi-fixed point where the solutions of renormalization group equations are 
concentrated at the electroweak scale. 
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1 Introduction 



A rapid development of experimental high-energy physics over the last decades of the 
XX century gave impetus to intensive investigations of various extensions of the Stan- 
dard Model. Its supersymmetric generalization known as the Minimal Supersymmetric 
Standard Model (MSSM) is one of the most popular extensions of the Standard Model. 
The Higgs sector of the MSSM includes two doublets of Higgs fields, Hi and H2. Upon a 
spontaneous breakdown of gauge symmetry, each of them develops a vacuum expectation 
value; we denote the corresponding vacuum expectation values by Vi and V2- The sum 
of the squares of the vacuum expectation values of the Higgs fields is v"^ = (246 GeV)^, 
the ratio of the expectation values being determined by the angle (3. By definition, 
f3 = arctan(f2/f 1). The value of tan/3 is not fixed experimentally. It is varied within a 
wide interval, from 1.3 — 1.8 to 50 — 60. Within supersymmetric (SUSY) models, the upper 
and lower limits on tan/3 arise under the assumption that perturbation theory is appli- 
cable up to the scale at which gauge coupling constants are unified, Mx = 3 ■ 10^^ GeV 
- that is, under the assumption that there is no Landau pole in solutions to relevant 
renormalization group equations. 

The spectrum of the Higgs sector of the MSSM contains four massive states. Two of 
them are CP-even, one is CP-odd, and one is charged. The presence of a light Higgs 
boson in the CP-even sector is an important distinguishing feature of SUSY models. Its 
mass is constrained from above as 



rrih < 



M| cos2 2/5 + A 



where Mz is the Z-boson mass {Mz ~ 91.2 GeV) and A stands for the contribution of 
loop corrections. The magnitude of these corrections is proportional to {rrit is the 
running mass of the t-quark) , depends logarithmically on the supersymmetry breakdown 
scale Ms, and is virtually independent of the choice of tan f3. An upper limit on the mass 
of the light CP-even Higgs boson within the MSSM grows with increasing tan (3 and, for 
tan/5 > 1, reaches 125 - 128 GeV in realistic SUSY models with Ms < 1000 GeV. 

At the same time it is known from ^ that, for tan/5 ^ 50 — 60, solutions to the 
renormalization group equations for the t-quark Yukawa coupling constant ht{t) are 
concentrated in the vicinity of the quasi-fixed point 
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6F{to) 
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with Qi being the gauge constants of the Standard Model group. The variable t is 
defined in the standard way: t = ln(M^/g^). Its value at the electroweak scale is 
to = 21n(Mx/MP°^'), where Mf°^' ^ 174.3 ± 5.1 GeV is the pole mass of the t-quark. 
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Along with the t-quark Yukawa couphng constant, solutions to the renormalization group 
equations for the corresponding trilinear coupling constant At characterising the inter- 
action of scalar fields and the combination DJlf = roq + m\j + m\ of the scalar particle 
masses also approach the infrared quasi-fixed point. The properties of solutions to the 
renormalization group equations within the MSSM and the spectrum of particles in the 
infrared quasi-fixed point regime at tan/3 ~ 1 were investigated in j2],|S]- 

A reduction of the number of independent parameters in the vicinity of the infrared 
quasi-fixed point considerably increased the predictive power of the theory. On the basis 
of the equation relating the Yukawa coupling constant for the t-quark with its mass at 
the electroweak scale 

mt[Mt ) = -= — fsm/i), (3) 



and the value calculated for the running mass of the t-quark within the MS'-scheme 
(mt(Mf°'^) = 165 ± 5 GeV), it was shown in 12-0 that, for a broad class of solutions 
satisfying the renormalization group equations within the MSSM and corresponding to 
the infrared quasi-fixed point regime, tan/5 takes values in the interval between 1.3 and 
1.8. These comparatively small values of tan/3 lead to much more stringent constraints on 
the mass of the lightest Higgs boson. A detailed theoretical analysis performed in 
revealed that, in the case being considered, its mass does not exceed 94 ± 5 GeV, which 
is 25 — 30 GeV below the absolute upper limit in the minimal SUSY model. It should 
be noted that the LEP II constraints on the mass of the lightest Higgs boson jOj are such 
that a considerable fraction of solutions approaching a quasi-fixed point at tan/? ~ 1 
have already been ruled out by existing experimental data. 

All the aforesaid furnishes a sufficient motivation for studying the Higgs sector in more 
complicated SUSY models, as well renormalization group equations and solutions to these 
equations therein. The present article is devoted to an analysis of coupling constant 
renormalization within a modified Next-to-Minimal SUSY Model (MNSSM), where the 
mass of the lightest Higgs boson can be as large as 120 — 130 GeV even at comparatively 
small values of tan/5 ~ 2. In addition, the spectrum of superpartners of observable 
particles and of Higgs bosons is studied in the vicinity of the quasi-fixed point of the 
renormalization group equations within the MNSSM. 



2 Modified NMSSM 

The Next-to-Minimal Supersymmetric standard Model (NMSSM) [7j-[n] is the simplest 
extension of the MSSM. Historically, the NMSSM arose as one of the possible solutions 
to the problem of the /i-term in supergravity (SUGRA) models [7]. Along with observ- 
able superfields, SUGRA theories contain a hidden sector that includes the dilaton and 
moduli fields {S and T^,, respectively), which are singlet in gauge interactions. The total 
superpotential in SUGRA models is usually represented as an expansion in the superfields 
of the observable sector; that is, 

W = 1^0(5, T„) + fi{S,Tj{HiH2) + ht{S,Tj{QH2)lf'R + ..., (4) 

where Wo{S,Tm) is the superpotential of the hidden sector. The expansion in (jH) pre- 
sumes that the parameter n appearing in front of the bilinear term {H1H2) must be about 
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the Planck mass, since this scale is the only dimensional parameter characterising the 
hidden sector of the theory. In this case, however, the Higgs bosons Hi and H2 acquire 
an enormous mass {m'jj_^ ~ /i^ ^ ^pi) and SU{2) ® f/(l) symmetry remains unbroken. 

In the NMSSM, an additional singlet superfield Y is introduced, while the term n{HiH2) 
is replaced by XY{HiH2) + (x/S)!''^. A spontaneous breakdown of gauge symmetry 
leads to the emergence of the vacuum expectation value {¥) = y/V^ of the field Y 
and to generation an effective //-term (/x = X{Y)). The resulting superpotential of the 
nonminimal SUSY model is invariant under discrete transformations of the group |HI. 
The Z3 symmetry of the superpotential of the observable sector naturally arises in string 
models, where all observable fields are massless in the limit of exact supersymmetry. 

Upon the introduction of the neutral field Y in the superpotential of the NMSSM, there 
arises the corresponding F-term in the potential of interaction of Higgs fields. As a 
result, an upper limit on the mass of the lightest Higgs boson becomes higher than that 
in the MSSM: 

mh< y yf 2 sin^ 2f3 + M| cos^ 2(3 + A . (5) 

In the tree approximation (A = 0), relation (0) was obtained in jH]. For A 0, the 
expressions for the upper limit in the MSSM and in the NMSSM coincide, after the 
substitution of Xy/\/2 for /i. The Higgs sector of the nonminimal SUSY model and one- 
loop corrections to it were studied in jTU],!!!]. In ^2]; the upper limit on the mass of the 
lightest Higgs boson within the NMSSM was contrasted against the analogous limits in 
the minimal standard and the minimal SUSY model. 

From relation it follows that the upper limit on rrih grows with increasing A(to)- It 
should be noted that only in the region of small values of tan f3 does it differ significantly 
from the analogous limit in the MSSM. As to the small tan jS scenario, it is realised in the 
case of sufficiently large values of htito). The growth of the Yukawa coupling constants 
at the electroweak scale is accompanied by an increase in /it(0) and A(0) at the Grand 
Unification scale; therefore, the upper limit on the mass of the lightest Higgs boson in 
the nonminimal SUSY model attains a maximum value in the limit of strong Yukawa 
coupling, in which case both /i^(0) and A^(0) are much greater than (7j^(0). 

Unfortunately, we were unable to obtain a self-consistent solution in the regime of strong 
Yukawa coupling within the NMSSM featuring the minimal set of fundamental parame- 
ters. Moreover, Z3 symmetry, which makes it possible to avoid the problem of the /x-term 
in the nonminimal SUSY model, leads to the emergence of three degenerate vacua in the 
theory. Immediately following the phase transition at the electroweak scale, the Universe 
is filled with three degenerate phases that must be separated by domain walls. However, 
the hypothesis of a domain structure of the vacuum is at odds with data from astro- 
physical observations. An attempt at destroying Z3 symmetry and the domain structure 
of the vacuum by introducing nonrenormalizable operators in the NMSSM Lagrangian 
leads to the appearance of quadratic divergences - that is, to the hierarchy problem ^13^ . 

In order to avoid the domain structure of the vacuum and to obtain a self-consistent 
solution in the regime of strong Yukawa coupling, it is necessary to modify the nonminimal 
SUSY model. The simplest way to modify the NMSSM is to introduce additional terms 
in the superpotential of the Higgs sector, fi{HiH2) and fi'Y^ 14 , that are not forbidden 
by gauge symmetry. The additional bilinear terms in the NMSSM superpotential destroy 
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^3 symmetry, and domain walls are not formed in such a theory. Upon the introduction 
of the parameter fi, it becomes possible to obtain the spectrum of SUSY particles in the 
modified model; for a specific choice of fi', the mass of the lightest Higgs boson reaches 
its upper limit, taking the largest value at x = 0. In analysing the modified NMSSM, 
it is therefore reasonable to set the coupling constant for the self-interaction of neutral 
superfields Y to zero. 

The MNSSM superpotential involves a large number of Yukawa coupling constants. At 
tan/3 ~ 1, they are all negligibly small, however, with the exception of the t-quark 
Yukawa coupling constant ht and the coupling constant A, which is responsible for the in- 
teraction of the superfield Y with the doublets Hi and H2- Thus, the total superpotential 
of the modified NMSSM can be represented in the form 

W^MNSSM = KH1H2) + l^'Y^ + XY{HiH2) + ht{QH2)UR • (6) 

Within SUGRA models, the terms in the superpotential ® that are bilinear in the super- 
fields can be generated owing to the term {Z{HiH2) + Z'Y"^ + h.c) in the Kahler potential 
or owing to the nonrenormalized interaction of fields from the observable and 
the hidden sector (this interaction may be due to nonperturbative effects) [THj.jTTj. 

Along with the parameters /i and /i', the masses of scalar fields mf, and the gaugino 
masses Mj are also generated upon a soft breakdown of supersymmetry. Moreover, a 
trilinear coupling constant Ai for the interaction of scalar fields is associated in the 
total Lagrangian of the theory with each Yukawa coupling constant, while a bilinear 
coupling constant B{B') is associated there with the parameter The hypothesis of 

universality of these constants at the scale Mx makes it possible to reduce their number 
to four: the scalar particle mass mg, the trilinear coupling constant A and the bilinear 
coupling constant Bq for the interaction of scalar fields, and the gaugino mass Mi/2- 

3 Analysis of the evolution of Yukawa couplings and 
determination of the quasi— fixed point 

The MNSSM parameters 

A, ht, fx, fi' , mo. A, Bo, M1/2 

specified at the Grand Unification scale evolve down to the electroweak scale or the 
scale of supersymmetry breakdown. Their renormalization is determined by the set of 
renormalization group equations, these equations for the coupling constants A, ht, Ai, mf, 
and Mi being coincident with the corresponding renormalization group equations within 
the NMSSM (see, for example, JT]) if one sets x = in them. The equations describing 
the evolution of /i, yu', B, and B' within the modified NMSSM were obtained in P^ j. 

Even in the one-loop approximation, the full system of renormalization group equations 
is nonlinear, so that it is hardly possible to solve it analytically. This set of equations 
can be broken down into two subsets. The first subset includes equations that describe 
the evolution of gauge and Yukawa coupling constants and of parameters /i and /i'. The 
second subset comprises equations for the parameters of a soft breakdown of supersym- 
metry. 
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In studying the evolution of the Yukawa couphng constants, it is convenient to introduce, 
instead of the constants ht, X, and Qi, the ratios 



^'^'~&Jt)' ^'^'~dJt)' ^'^'~dJt)' 

where Yx{t) = A^(t)/(47r)^. The region of admissible values of the Yukawa coupling 
constants at the electroweak scale is bounded by the quasi-fixed (or Hill) line. Beyond 
this region, solutions to the renormalization group equations for Yi(t) develop a Landau 
pole below the scale Mx, so that perturbation theory becomes inapplicable for ~ ^x- 
The results of our numerical calculations are presented in Fig. 1, whence one can see that, 
in the regime of strong Yukawa coupling, all solutions for Yi{t) are attracted to the Hill 
line, which intersects the pt axis at the point whose coordinates {p\{to), ptitoj) = (0, 0.087) 
correspond to the quasi-fixed point in the MSSM. 

In analysing the results of the numerical calculations (see Fig. 1), attention is captured by 
a pronounced nonuniformity in the distribution of solutions to the renormalization group 
equations along the quasi-fixed line. The main reason behind this is that, in the regime 
of strong Yukawa coupling, solutions are attracted not only to the quasi-fixed but also to 
the infrared fixed (or invariant) line. The latter connects two fixed points. One of them 
is the stable infrared fixed point for the set of renormalization group equations within 
the MNSSM {pt = 7/18, pa = 0,pi = 0,p2 = 0) As the invariant line approaches 

this point, px ~ (pt ~ 7/18)^^/^^. The other fixed point {p\/ Pt) = 1 corresponds to large 
values of the Yukawa coupling constants, Yt,Y\ ^ cij, in which case the gauge coupling 
constants can be disregarded [THj. In the limiting case of p\, pt ^ 1, the asymptotic 
behaviour of the curve being studied is given by 

PA = - - - ^^P.. (7) 

The infrared fixed lines and their properties in the minimal standard and the minimal 
supersymmetric model were studied in detail elsewhere 



With increasing initial values of the Yukawa coupling constants Yt{Q) and Y\{Q) at the 
Grand Unification scale, the region where solutions are concentrated at the electroweak 
scale shrinks abruptly and all solutions to the renormalization group equations within the 
MNSSM are focused near the point of intersection of the invariant and the quasi-fixed 
line: 

p?^^(to) = 0.803, pfP(to) = 0.224. (8) 

This point can be considered as the quasi-fixed point for the set of renormalization group 
equations for the modified NMSSM 



Among subsidiary constraints that are frequently imposed in studying supersymmetric 
models, we would like to mention the unification of the Yukawa coupling constants for 
the 6-quark and the r-lepton at the scale Mx'-, this usually occurs in minimal schemes for 
unifying gauge interactions - for example, in those that are based on the SU{b), the E^, 
or the 5*0(10) group. The unification of /if, and h^- within the MNSSM is realised only 
in the case where the constants satisfy a specific relation between Y^ and Yx- Integrating 
the renormalization group equations and substituting Rhrito) = m;,(to)/'n^r(^o) = 1-61, 
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which corresponds to mrito) = 1.78 GeV and mfe(to) = 2.86 GeV, we obtain 
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3.67 



YM 



(9) 

The results obtained here indicate that b — t unification is possible under the condition 
that Yt{0) ^ Ytito), which is realised only in the regime of strong Yukawa coupling. By 
varying the running mass mh{mh) of the 6-quark from 4.1 to 4.4 GeV, we found that the 
equality of the Yukawa coupling constants at the Grand Unification scale can be achieved 
only at tan f3 < 2. 

The possibility of unifying the Yukawa coupling constants within the NMSSM was com- 
prehensively studied in The condition 1^6(0) = 5^r(0) imposes stringent con- 
straints on the parameter space of the model being studied. Since hb and hr are small 
in magnitude at tan/? ~ 1, they can be generated, however, owing to nonrenormalizable 
operators upon a spontaneous breakdown of symmetry at the Grand Unification scale. In 
this case, hb and h^. may be different. In studying the spectrum of superpartners below, 
we will not therefore assume that RbriO) = 1. 



4 Renormalization of the soft SUSY breaking pa- 
rameters 

If the evolution of gauge and Yukawa coupling constants is known, the remaining subset 
of renormalization group equations within the MNSSM can be treated as a set of linear 
differential equations for the parameters of a soft breakdown of supersymmetry. For 
universal boundary conditions, a general solution for the trilinear coupling constants 
Ai{t) and for the masses of scalar fields m'^it) has the form 

A(t)=e,(t)A + /,(t)Mi/2, 
Tfiit) = ai{t)ml + h{t)Ml/^ + Ci{t)AMi/2 + di{t)A^ . 

The functions ej(t), aiit), bi(t), Ci{t), and di(t), which determine the evolution of 

Ai{t) and m^{t), remain unknown. The results of our numerical calculations reveal that 
these functions greatly depend on the choice of Yukawa coupling constants at the scale 
Mx. 

In analysing the behaviour of solutions to the renormalization group equations in the 
regime of strong Yukawa coupling, it is more convenient to consider, instead of the squares 
of the scalar particle masses, their linear combinations 

ml{t)=ml{t)+ml{t)+ml{t), 

ml{t) = m\{t) + ml{t) + mlit) . ^ ' 

For the universal boundary conditions, solutions to the differential equations for 93Tf(t) 
can be represented in the same form as the solutions for m^(t) (see (|T0|) ): that is 

Tfiit) = 3diit)ml + h{t)Ml/^ + Ci{t)AMi/2 + di{t)A^. (12) 
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Since the homogeneous equations for Ai{t) and have the same form, the functions 

di{t) and ej(t) coincide. 

With increasing 1^(0), the functions ej(to), Cj(to), and (ij(to) decrease and tend to zero in 
the hmit Yi{0) oo. Concurrently, At(t), Ax{t), DJl^it), and become independent 

of A and rriQ, while relations pO|) and p2|) are significantly simplified. This behaviour of 
the solutions in question implies that, as the solutions to the renormalization group equa- 
tions for the Yukawa coupling constants approach quasi-fixed points, the corresponding 
solutions for Ai(t) and also approach the quasi-fixed points whose coordinates are 



pSf (to)^1.77, pg^f (to)^6.09, 

(13) 



pfjito) ^ -0.42, pg|P(to) ^ -2.^ 

where pA.it) = Ai{t)/Mi/2 and p<xs{^{t) = DJlf/M^^^. At the same time, the functions ai{t) 
approach some constants independent of t and Yi{0): 

a,(t)^l/7, ai(t)^4/7, a,it) ^ -5/7, 

a„(t)^l/7, a,(t)->4/7. ^ ' 

In the case of nonuniversal boundary conditions at Yt{0) ~ ^a(O), the required solution 
to the renormalization group equations for Ai{t) and OJtf (t) can be represented as 



(15) 



where and /?i are constants of integration that can be expressed in terms of At{0), 
Aa(0), 9712(0), and aJli(O); ei(t) = Yt{t)/Yt{0); Ai = 1 and As = 3/7. The functions Vij{t) 
are weakly dependent on the Yukawa coupling constants at the scale Mx, and Vij{0) = 1. 
In equations (fT3j) . we have omitted terms proportional to M1/2, ^1/2; ^j(0)Mi/2, and 
A(0)A,(0). 

With increasing Vt(0) ^ yx(0), the dependence of v4j(to) and 9Jl^(to) on ai and Pi quickly 
becomes weaker. The results of our numerical analysis that are displayed in Fig. 2 
indicate that, for hf{0) = X^i^O) = 20 and boundary conditions uniformly distributed in 
the {At, A\) and the (971^, DJl\) plane, the solutions to the renormalization group equations 
for the parameters of the soft SUSY breaking in the vicinity of the quasi-fixed point 
are concentrated near some straight lines. The equations of these straight lines can be 
obtained by setting Aa(0) = -3A(0) and ajt|(0) = -39712 (0) ^^-^^^ = 0) 

at the Grand Unification scale. As a result, we find that, at the electroweak scale, the 
parameters of a soft breakdown of supersymmetry satisfy the relations 



A + 0.137(0.147)Aa = I.7OM1/2, 

(16) 

/2 



mf + 0.137(0.147)971^ = 5.76Ml^. 



The equation for 971^ has been obtained for all Ai{0) set to zero. In relations (fT^ . the 
coefficients obtained by fitting the results of our numerical calculations (see Fig. 2) are 
indicated parenthetically. As the Yukawa coupling constants approach quasi-fixed points, 
the two sets of coefficients in ()16|) approach fast each other and, at Yi{0) ~ 1, become 
virtually coincident. 
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5 Spectrum of SUSY particles and Higgs bosons 



Let us now proceed to study of the spectrum of the superpartners of observable particles 
and Higgs bosons in the vicinity of the quasi-fixed point within the MNSSM. The Yukawa 
coupling constants ht and A are determined here by relations (jH))- The value of tan/3 
can be calculated by formula (jHl). In the regime of the infrared quasi-fixed point at 
mt{Mf"^'') = 165 GeV the result is tan/5 = 1.88. 

The remaining fundamental parameters of the MNSSM must be chosen in such a way 
that a spontaneous breakdown of SU{2)^U{1) gauge symmetry occur at the electroweak 
scale. The position of the physical minimum of the potential representing the interaction 
of Higgs fields is determined by solving the set of nonlinear algebraic equations 

dV{vi,V2,y) ^ Q dV{vuV2,y) ^ ^ dV{vi,v2,y) ^ ^ ^^^^ 
dvi ' dv2 ' dy ' 

where V{vi,V2,y) is the effective potential of interaction of Higgs fields within the 
MNSSM nil. 

Since the vacuum expectation value v and tan (3 are known, the set of equations (fT7j) can 
be used to determine the parameters fi and Bq and to compute the vacuum expectation 
value {Y). Instead of /x, it is convenient to introduce here ^cs = + Ay/-\/2. The sign 
of /ieflf is not fixed in solving the set of equations (fTTj) : it must be considered as a free 
parameter of the theory. The results obtained in this way for the vacuum expectation 
value y, the parameters fi^s and Bq, and the spectrum of particles within the modified 
NMSSM depend on the choice of A, mo, M1/2, and fi'. 

It is of particular interest to analyse the spectrum of particles in that region of the 
parameter space of the MNSSM where the mass of the lightest Higgs boson is close to 
its theoretical upper limit, since the remaining part of the parameter space is virtually 
ruled out by the existing experimental data. For each individual set of the parameters 
A, mo, and M1/2, the mass of the lightest Higgs boson reaches the upper bound on itself 
at a specific choice of fi'. It is precisely at these values of the parameter fi' that we 
have calculated the particle spectrum presented in Tables 1 and 2. On the basis of our 
numerical results given there, one can judge the effect of the fundamental constants A, 
mo, and M1/2 on the spectrum of the superpartners of the t-quark (m^-^ 2)' gluino 
(M3), of the neutralino (^«xi,...,5)' chargino (m-±^), and of the Higgs bosons {m^, 

f^H, fns, ^Ai^2)- each set of the aforementioned parameters, we quote the mass 
of the hghtest Higgs boson according to the calculations in the one- and the two-loop 
approximation, along with the corresponding values of /lefr, Bo, y, and fi'. As can be seen 
from the data displayed in Tables 1 and 2, the qualitative pattern of the spectrum within 
the MNSSM undergoes no changes in response to variations of the parameters A and mo 
within reasonable limits. 

The CP-even Higgs boson (m^), which corresponds to the neutral field Y is the heavi- 
est particle in the spectrum of the modified NMSSM, while the neutralino (jt^xs) 
heaviest fermion there, the main contribution to its wave function coming from the super- 
partner of the field Y. With increasing mg the masses of the squarks, the Higgs bosons, 
and the heavy chargino and neutralino grow, whereas the spectrum of extremely light 
particles remains unchanged. Since the dependence of the parameters of a soft break- 
down of supersymmetry on A disappears at the electroweak scale in the regime of strong 
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Yukawa coupling, the parameters fies, B, and /x', together with the spectrum of the su- 
perpartners of observable particles and the mass of the lightest Higgs boson, undergo 
only slight changes in response to a variation of the trilinear coupling constant for the 
interaction of scalar fields from —M1/2 to Mi/2- Despite this, the A dependence of masses 
of one of the CP-even (ms) and two CP-odd ("^^1.2) Higgs bosons survives. It is due 
primarily to the fact that the bilinear coupling constant B' for the interaction of neutral 
scalar fields is proportional to A. It should be noted in addition that, for a specific choice 
of fundamental parameters (in particular, of the parameter A), the mass of the lightest 
CP-odd Higgs boson may prove to be about 100 GeV or less. However, this Higgs boson 
takes virtually no part in electroweak interactions, since the main contribution to its 
wave function comes from the CP-odd component of the field Y. Therefore, attempts at 
experimentally detecting it run into problems. 

Loop corrections play an important role in calculating the mass of the lightest Higgs 
boson. Their inclusion results in that the mass of the lightest Higgs boson proves to be 
greater for fi^s < than for /igfr > 0. This is because rrih grows as the mixing in he sector 
of the superpartners of the t-quark {tn and ti) becomes stronger. The point is that the 
mixing of in and is determined by the quantity Xt = At + fies/ tan/? and is therefore 
greater in magnitude for fies < since At < 0. It should also be noted that the inclusion 
of two-loop corrections leads to a reduction of m/j by approximately 10 GeV. The mass 
of the lightest Higgs boson depends only slightly on A and mo, because the squark masses 
depend slightly on the corresponding fundamental parameters (see Tables 1 and 2). The 
value of rrih is determined primarily by the supersymmetry breaking scale Ms - that is, 
by the quantity M3. From our numerical results quoted in Tables 1 and 2, one can see 
that, at mt{Mf"^^) = 165 GeV and M3 < 2 TeV, the mass of the lightest Higgs boson 
does not exceed 127 GeV. 



6 Conclusions 

We have studied coupling constant renormalization and the spectrum of particles within 
the modified Next-to-Minimal Supersymmetric Standard Model (MNSSM). We have 
shown that, in the regime of strong Yukawa coupling, solutions to the renormalization 
group equations for Yi(t) are attracted to the Hill line and that, under specific conditions 
b — T unification is realised at the scale Mx- In the limit ^i(O) — >■ 00, all solutions for the 
Yukawa coupling constants are concentrated in the vicinity of the quasi-fixed point that 
is formed in the space of the Yukawa coupling constants as the result of intersection of 
the invariant and the Hill line. 

As the Yukawa coupling constants approach the quasi-fixed point, the corresponding 
trilinear coupling constants and combinations ()11|) of the scalar particle masses cease to 
depend on the boundary conditions at the scale Mx- In the case of nonuniversal boundary 
conditions, Ai(t) and ^j(t) are attracted to straight lines in the space spanned by the 
parameters of a soft breakdown of supersymmetry and, with increasing ^^(O), approach 
the quasi-fixed points, moving along these straight lines. 

We have analysed the spectrum of particles in the infrared quasi-fixed point regime of 
the MNSSM. The CP-even Higgs boson, which corresponds to the neutral field Y, is the 
heaviest particle in this spectrum. At reasonable values of the parameters of the model 
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being studied, the gluinos, the squarks, and the heavy Higgs bosons are much heavier 
than the hghtest Higgs boson and than the hghtest chargino and the hghtest neutrahno as 
well. This is not so only for one of the CP-odd Higgs bosons whose mass changes within 
a wide range in response to variations in the fundamental parameters of the MNSSM. In 
the vicinity of the quasi-fixed point at mt(Mf°''^) = 165 GeV, the mass of the lightest 
Higgs boson does not exceed 127 GeV. 
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Table 1. Particle spectrum in the vicinity of the quasi-fixed point in the MNSSM 
at mt(Mf°'^) = 165 GeV, tan/3 fa 1.883, and /leS > depending upon the choice of 
fundamental parameters A, mo, and M1/2 (all parameters and masses are given in GeV). 





ml 





Ml,2 














A 








-Ml/2 


O.5M1/2 








Ml/2 


-392.8 


-392.8 


-392.8 


-392.8 


-785.5 


-196.4 


mt{to) 


165 


165 


165 


165 


165 


165 


tan/3 


1.883 


1.883 


1.883 


1.883 


1.883 


1.883 


f^eS 


728.6 


841.7 


726.8 


730.1 


1361.2 


380.4 


Bo 


-1629.1 


-1935.4 


-1260.0 


-1813.2 


-3064.4 


-861.8 


y 


-0.00037 


-0.00021 


-0.00043 


-0.00035 


-0.00006 


-0.00233 




-1899.8 


-2176.7 


-1905.9 


-1898.3 


-3544.6 


-993.1 


mh(to) 


125.0 


125.1 


125.0 


125.0 


134.9 


114.8 


(1-loop) 














mh(to) 


118.4 


118.5 


118.4 


118.4 


123.2 


111.9 


(2-loop) 














M3(l TeV) 


1000 


1000 


1000 


1000 


2000 


500 


mi,(l TeV) 


840.6 


889.7 


841.1 


840.3 


1652.0 


447.4 


m,-,(l TeV) 


695.1 


713.6 


696.6 


694.3 


1366.2 


371.6 


mffil TeV) 


898.5 


1080.5 


895.4 


900.3 


1691.0 


468.8 


m5(l TeV) 


2623.4 


3034.3 


2452.2 


2706.0 


4901.7 


1378.0 


mA,{l TeV) 


953.9 


1113.8 


1245.7 


925.2 


1722.6 


538.2 


mA,{l TeV) 


704.3 


762.7 


872.0 


318.2 


1366.2 


302.2 




164.6 


164.4 


164.6 


1 a A a 

164.6 


326.9 


OA 

84.3 


m.^,{to) 


327.8 


327.6 


327.8 


327.8 


649.4 


170.1 


m^3(l TeV) 


755.1 


870.8 


753.3 


756.7 


1404.2 


400.9 


1^x4(1 TeV)| 


755.9 


872.6 


755.1 


758.4 


1405.0 


404.3 


m^,(lTeV)| 


1931.8 


2212.3 


1938 


1930.3 


3599.0 


1015.4 


m-±(to) 


327.8 


327.6 


327.8 


327.8 


649.4 


169.9 


m-±(lTeV) 


757.0 


872.6 


755.2 


758.5 


1405.2 


404.5 
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Table 2. Particle spectrum in the vicinity of the quasi-fixed point in the MNSSM 
at mt(Mf°'^) = 165 GeV, tan/3 fa 1.883, and /leS < depending upon the choice of 
fundamental parameters A, mo, and M1/2 (all parameters and masses are given in GeV). 



ml 





^1/2 














A 








-Ml/2 


Ml/2 








Ml/2 


-392.8 


-392.8 


-392.8 


-392.8 


-785.5 


-196.4 


mt{to) 


165 


165 


165 


165 


165 


165 


tan P 


1.883 


1.883 


1.883 


1.883 


1.883 


1.883 




-727.8 


-840.9 


-726.0 


-731.2 


-1360.7 


-378.9 


Bo 


1008 


1320.3 


1366.7 


647.9 


2050.4 


495.8 


y 


-0.00149 


-0.001 


-0.00128 


-0.00177 


-0.00020 


-0.0112 


/^'(to) 


1671.5 


1950.6 


1656.8 


1690.3 


3172.7 


857.8 


mh(to) 


134.1 


134.9 


134.0 


134.2 


143.1 


124.1 


(1-loop) 














mh(to) 


124.4 


124.8 


124.3 


124.5 


127.2 


119.6 


(2-loop) 














M3(l TeV) 


1000 


1000 


1000 


1000 


2000 


500 


m,-^(lTeV) 


890.2 


935.6 


890.5 


889.8 


1682.8 


507.9 


^t.(lTeV) 


630.3 


652.2 


632.2 


628.0 


1328.1 


283.5 


mj^(l TeV) 


896.2 


1078.5 


893.5 


899.3 


1689.9 


464.4 


m5(l TeV) 


2147.4 


2565.9 


2309.2 


1972.3 


4126.5 


1097.7 


m^,(lTeV) 


1123.2 


1219.3 


931.0 


1437.9 


1984.8 


623.1 


mA^il TeV) 


857.6 


1017.8 


545.0 


886.9 


1657.5 


412.8 


m^,{to) 


160.0 


160.5 


160.0 


160.0 


324.4 


74.9 


m^^{to) 


311.1 


313.7 


311.0 


311.2 


639.9 


141.4 


\m^,{l TeV)| 


753.7 


896.6 


751.9 


757.2 


1403.4 


398.5 


m^^{l TeV) 


764.7 


878.1 


763.0 


768.1 


1410.0 


416.7 


m^,(l TeV) 


1700.7 


1983.2 


1685.8 


1719.6 


3221.8 


879.1 


m-±(to) 


310.7 


313.4 


310.7 


310.8 


639.8 


139.4 


m-±(lTeV) 


763.3 


877.0 


761.6 


766.7 


1409.1 


AU.5 
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Figure captions 



Fig. 1. Boundary conditions for the renormalization group equations of the MNSSM 
at the scale q = Mx uniformly distributed in the {pt, p\) plane in a square 
2 < /i^(0),A^(0) < 10 - Fig. la, and the corresponding values of the Yukawa cou- 
plings at the electroweak scale - Fig. lb. The thick and thin curves in Fig. lb represent, 
respectively, the invariant and the Hill line. The dashed straight line in Fig. lb is a fit of 
the values {pt{to), px{to)) for 20 < h^{0),X'^{0) < 100. 

Fig. 2. Boundary conditions for the renormalization group equations of the MNSSM 
at the Grand Unification scale {t = 0) at hf{0) = A^(0) = 20 uniformly distributed in 
the {At/Mi/2, A\/Mi/2) plane - Fig. 2a, and the corresponding values of the trilinear 
couplings at the electroweak scale {t = to) - Fig. 2b. The straight line in Fig. 2b is a fit 
of the values {At{to) , Ax{to)) . 

Fig. 3. Boundary conditions for the renormalization group equations of the MNSSM 
at the Grand Unification scale {t = 0) at /if(0) = X^iO) = 20 and ^^(0) = ^^(0) = 
uniformly distributed in the {^M^/My^, dJll/My^) plane - Fig. 3a, and the corresponding 
values of the trilinear couplings at the electroweak scale {t — to) - Fig. 3b. The straight 
line in Fig. 3b is a fit of the values {Tl^{to),Tll{to)). 
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Fig.la. 



Fig.lb. 



